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Abstract 

The theory of spin superradiance, developed earher for nuclear magnets, is general- 
ized to a wider class of spin systems, such as granular magnets and molecular magnets. 
The latter may possess strong single-site magnetocrystalline anisotropy, whose role in 
nonlinear spin dynamics is analysed. Transient as well as pulsing superradiant regimes 
are described. These coherent regimes may be employed in the operation of spin masers. 
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1 Introduction 



Among different masers, generating radiation at microwave and radio frequencies [1], tliere 
is a separate class of spin masers, in wliicli the radiation process is due to moving spins. Tlie 
role of a resonator for such masers is played by a resonant electric circuit coupled with a spin 
system. The peculiarity of spin masers is that their magnetodipole radiation is rather weak 
and practically does not propagate into free space but mainly is taken up by a resonant coil 
surrounding the spin sample [2]. Nevertheless, because of a deep physical similarity with 
other types of masers, radiating spin systems are called spin masers [3]. 

Resonant optical systems can display the effect of superradiance, which is a coherent spon- 
taneous emission, when the radiation intensity is approximately proportional to the number 
of radiators squared [4,5]. Spin systems may also exhibit this phenomenon [6] which, because 
of its direct analogy with optical atomic superradiance, is termed spin superradiance. The 
latter, in the same way as the former, can be of two major types, transient and pulsing. The 
transient superradiance occurs as a single sharp supcrradiant burst that may be accompa- 
nied by several small quickly diminishing oscillations. This regime happens when the system 
is prepared in an inverted strongly nonequilibrium state, after which it is not influenced 
by any additional external fields. The pulsing superradiance corresponds to a long train of 
supcrradiant bursts, which can be realized if the system is subject to a permanent pump- 
ing mechanism. Both types of superradiance were experimentally observed for several spin 
systems, including the transient [7-10] and pulsing [11-13] spin superradiance. The theory 
of spin superradiance was developed [14-20], being based on microscopic Hamiltonians and 
allowing for the first correct description of purely self- organized superradiance. Numerical 
simulations of spin superradiance, being a kind of computer experiments, were also realized 
[21-27]. 

These studies of spin superradiance were based on the Hamiltonian describing an ensem- 
ble of polarized nuclear spins interacting through dipolar forces. Such a matter, for short, 
can be called a nuclear magnet. The influence of hyperfine forces on nuclear spin superra- 
diance was also considered [28-30]. This attention to nuclear magnets has been due to the 
fact that the experiments with spin superradiance [7-13] had been accomplished for nuclear 
spins. In principle, this effect could be realized for electronic spins as well, with a similar 
theoretical description of electron spin superradiance [28,31,32]. 

In the present paper, the theory of spin superradiance is generalized to a wider class of 
materials, such as granular magnets and molecular magnets, which can possess higher spins 
and strong single-site magnetic anisotropy. The theory is also improved by taking a more 
accurate account of retardation effects. 

2 Materials Characteristics 

The experiments on spin superradiance till now have been accomplished with nuclear mag- 
nets. The first observation of transient spin superradiance [7,8] and its confirmation [9] were 
done for propanediol C3H8O2. An active substance here is the ensemble of proton spins with 
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the density ~ 4 x 10^^ cm~^. The spins were polarized, by means of dynamic nuclear po- 
larization, parallel to an external magnetic field So ~ 1 T, which corresponds to the Zeeman 
frequency a;o ~ 10^ Hz. The material was kept at low temperature T ~ 0.1 K, which resulted 
in strong suppression of the nuclear spin-lattice relaxation. The corresponding longitudinal 
relaxation time was Ti ~ 10"'' s. The transverse dephasing time, due to dipolar interactions, 
was T2 ~ 10~^ s. The sample was coupled to a resonant electric circuit with a quality factor 
Q ~ 100 and a ringing time r ~ 10~^ s. 

In other experiments on transient spin superradiance [10] butanol C4H9OH and ammo- 
nia NH3 were used. These materials are rich with protons of density pn ~ 10^^ cm~^. An 
external magnetic field i?o ~ 1 T defined the Zeeman frequency Uq ~ 10^ Hz. The experi- 
ments were carried out at low temperature, resulting in the long spin-lattice relaxation time 
Ti ~ 4 X 10^ s up to 10^ s. The transverse dephasing time was T2 ~ 10~^ s. The quahty 
factor of the resonant electric circuit was Q ~ 30, and the ringing time r ~ 5 x 10~^ s. 

The experiments on pulsing spin superradiance [11-13] employed the ruby crystal AI2O3. 
The active nuclei here are ^^Al, with spin 7 = 5/2 and density p^i ~ 4 x 10^^ cm^^. The 
crystal was oriented in an external magnetic field Bq r-^ IT so that a fully resolved structure 
of its five Am = ±1 transitions could be observed. Then, if a resonant circuit is tuned to 
a selected transition line, ^^Al spins form an effective two-level system. In experiments, the 
circuit was tuned to the central { — |, |} hne, with a transition frequency cjq ~ 10^ Hz. At 
low temperatures T ~ 1 K, the spin-lattice relaxation time was Ti ~ 10^ s. The transverse 
relaxation time was T2 ~ 10~^ s. The quality factor of the resonant circuit was Q ~ 100, with 
the ringing time r ~ 10~^ s. The inversion of spin polarization was permanently supported 
by means of dynamic nuclear polarization with the pumping rate — 10 s. 

Wishing to extend the possibility of realizing spin superradiance for other types of spin 
systems, the first such materials that come to mind are granular magnets. These are com- 
posed of magnetic nanoparticles of diameters between 10 A to 10^ A. Each nanoparticle is a 
superparamegnetic cluster of an effective spin S that can be sufficiently high. There exists a 
large variety of magnetic nanoparticles [33,34] which can be formed by simple metals, such 
as Ni, Fe, Co, and Hg or their oxides, as NiO and Fe203. Many magnetic nanoparticles 
are made of different alloys, such as NiFe204, Nd2Fei4B, Pr2Fei4B, Tb2Fei4B, Dy2Fei4B, 
Pr2Coi4B, SmiFciiTii, SmiFeioV2, Sm2Fei7N23, Sm2Fei7C22, Sm2Coi7, and SmCos. An 
ensemble of magnetic nanoparticles could be polarized in an external magnetic field, after 
which, inverting the latter, one would get an inverted nonequilibrium system. The following 
process should be similar to that developing in an ensemble of nuclear spins. The main disad- 
vantage of granular magnets is that composing them nanoparticles vary in size and shape. It 
is practically hardly possible to make a system of nanoparticles being almost identical. And 
an essential variation of the properties of radiating objects leads to a large inhomogeneous 
broadening, which hinders the feasibility of achieving a good level of coherence. 

Another class of composite objects, possessing nonzero spin, are magnetic molecules 
[35]. These molecules can form crystalline materials where all magnetic clusters are well 
defined with the same shape, size, and orientation, because of which the inhomogeneous 
broadening has to be very low. Such materials, composed of magnetic molecules are termed 
molecular magnets. There are many different magnetic molecules [36,37] having in their 
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ground state nonzero total spins S. For example, the dodecanuclear manganese cluster 
[Mni20i2(CH3COO)i6(H20)4]-2CH3COOH-4H20 is a molecule with spin S = 10. It has a 
rather strong single-site anisotropy characterized by a parameter D ^ 0.7 K, which makes the 
anisotropy barrier DS^ ~ 70 K. At low temperatures, lower than the blocking temperature 
Tb ~ 3 K, the magnetization of a molecular crystal is preserved during the relaxation time 
~ 10^ s. The size of a molecule is about 10 A and the distance between neighbouring 
molecules in the molecular crystal is about 14 A. The interaction between molecules is 
through dipolar forces, with an energy 0.1 K or 10^° Hz. More information on the 
properties of this molecule can be found in Refs. [38-48]. 

The octanuclear iron molecular cluster [Fe802(OH)i2(tacn)6]^''', where "tacn" stands for 
the organic ligand triazacyclononanc, also has a high spin S = 10. Its magnetic anisotropy 
is D ~ 0.3 K, hence the anisotropy barrier is DS^ ~ 30 K. The blocking temperature is 
Tb ~ 1 K, below which the relaxation time for a molecular crystal is about 10^ s. The 
size of a molecule and the intermolecular distance in a crystal are close to those for the 
molecule Mni2 mentioned above. Other characteristics of the molecule can be found in Refs. 
[36,37,49-53]. 

The molecule [Cr(CNMnL)6] (0104)9, where L is a neutral pentadentate ligand, has in 
its ground state the total spin S = 27/2 (see [37]). The molecule [(PhSi02)6Cu6(02SiPh)6] 

possesses the spin S = 3 (sec [54]). And the molecule K6[Vi^As6042(H20)]-8H20 has the 
spin S = 1/2 and no magnetic anisotropy [55,56]. There is a number of other molecules [37] 
with nonzero total ground-state spin. Some molecules may have zero spin in their ground 
state but a finite spin in excited states [57-59]. 

The relaxation of the total magnetization of a crystal, consisting of many magnetic 
molecules, occurs because of the axial degeneracy of spin direction in each molecule. At 
zero external magnetic field, the spin of a single molecule can be directed either up or down, 
so that the equilibrium state of an ensemble of these molecules corresponds to zero total 
magnetization. At temperatures higher than the blocking temperature Tg, the relaxation is 
rather fast and is due to thermal fluctuations. At low temperatures below Tb, the relaxation 
is very slow, so that the polarization of a molecular magnet can be blocked for months. At 
such low temperatures, the relaxation is characterized as quantum spin tunneling between 
the degenerate states. To understand the low-temperature behaviour of magnetic molecules, 
it is necessary to take into account both the internal effects, caused by atoms composing each 
molecule [60-62], as well as external interactions between the molecules forming a crystal 
[63-65]. Taking account of dipole interactions between molecules is crucial for correctly 
describing the relaxation in a molecular magnet [42,51,52,56,66]. 

In this way, there exists a large variety of different objects, such as nuclei, granules, and 
molecules, each of which can be considered as an entity, like a particle with an effective spin 
that can vary in a wide diapason. The main interactions between such magnetic particles, 
forming a solid, are dipole interactions. A sample, composed of these particles, can be 
polarized and at low temperatures the polarization can be preserved for a very long time. A 
specific feature of magnetic particles with high spin is the presence of the single-site magnetic 
anisotropy, which one has to take into account when considering collective spin dynamics. 
The problem, to be addressed in the following sections, is how coherent spin radiation can 
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arise is such materials made of objects with high spins and magnetic anisotropy. A special 
attention will be paid to the possibility of realizing spin superradiance that is a self-organized 
process developing without imposing on the spin system an initial transverse coherence. 

3 Types of Coherence 

Before passing to the development of a generalized theory of spin superradiance, a few words 
are to be said for concretizing the term "coherence" that will be repeatedly used in what fol- 
lows. Generally, one uses this term in two different meanings. One widespread usage implies 
under a coherent state of a many-particle system just a pure quantum state characterized by 
the same wave function for all particles. Then an incoherent state is a mixed state described 
by a density matrix [67]. In the theory of nuclear magnetic resonance [68,69], coherence 
usually means the existence of transverse magnetization and, more generally, the existence 
of non-diagonal matrix elements. One also tells [68,69] that the existence of the transverse 
magnetization means phase coherence, as opposed to amplitude coherence associated with 
a nonzero longitudinal polarization. To formalize these definitions for a system of N spins, 
let us introduce the notation 

1 N 1 ^ 

AT ^ * AT ^ 

-'^ i=l -'^ 1=1 

where is the ^-component of the spin operator and Sf are the raising and, respectively, 
lowering spin operators. Denote the statistical averaging by the angle brackets < ... >. 
Then a nonzero < S*^ means longitudinal coherence, or diagonal coherence, or amplitude 
coherence. While a nonzero < >^ signifies the existence of transverse coherence, or 
nondiagonal coherence, or phase coherence. 

Another possibility could be to tell that a nonzero < >^ corresponds to state 
coherence. When the whole system is in pure state, then < >= ±S, which can be 
named pure coherent state, to distinguish it from the case of partially coherent state, when 
< I < 5"^ > I < 5". This definition of state coherence is in agreement with the quantum- 
mechanical understanding of a coherent state as of a pure state. It also agrees with the 
definition of coherent states in quantum field theory as of eigenvalues of field operators [70]. 
In many applications, the operator characterizes a population difference. Hence, nonzero 
< S*^ > may be associated with population coherence. Being related to the determination of 
eigenfunctions of operators, the state coherence may also be called quantum coherence. 

Since the raising and lowering operators describe transitions between quantum states, 
the existence of nonzero < >7^ can be ascribed to transition coherence. This type of 
coherence is closely connected with the existence of coherent radiation, because of which it 
may be termed radiation coherence. Also, this recalls the classical understanding of coherence 
as of synchronous motion of several objects, which makes it sometimes possible to use the 
name of classical coherence. Synchronous motion is often termed as the motion in phase. 
That is why the term of phase coherence is appropriate here. 

These types of coherence are not strictly correlated with each other. A spin system may 
possess one of them or both, or neither. But often they are complimentary to each other, as 
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it happens in the process of superradiance. Then, at the initial time, the spin system has to 
be prepared in a well polarized state, thus, displaying state coherence. In addition, this state 
must be strongly nonequilibrium. Transition coherence, at the initial time, should be absent. 
It has to develop in a self- organized way owing to mutual spin correlations, which can be 
realized trough a feedback field. Usually, the maximal transition coherence in superradiance 
develops at the moment when the state coherence is minimal. 



4 Spin Hamiltonian 



Consider an ensemble of N spins Sj enumerated by the index i — 1,2, ... ,N. These spins 
can correspond either to nuclei, or to granules, or to molecules. The Hamiltonian of the 
system, 

H^^Hi + lT^Hij, (1) 

contains the terms Hi, related to individual spins, and the terms Hij describing spin inter- 
actions. An individual term 

Hi = -//oB • S, - D{Stf (2) 

includes the Zeeman energy and a part characterizing magnetic anizotropy. In the case 
of nuclei, //q = gs^Jl'N — ^75, where gs is the Lande factor; fxj^, nuclear magneton; 75, 
gyromagnetic ratio. For nuclei, /xq can be positive as well as negative. When S is an electron 
spin, //Q = —Qs^^B = —fi'ys, where fiB is the Bohr magneton. Then Hq is negative. The 
anisotropy term, with an anisotropy constant D, is nontrivial only for spins higher than 1/2. 
Positive D > implies an easy- axis anisotropy, while D < means an easy-plane anisotropy. 
The pair terms in the Hamiltonian (1) describe dipolar interactions 

Hij^Y^CffsrS^ (3) 

a/3 



through the dipolar tensor 

< - # (^"/^ - 3<<) , (4) 



in which a,(5 — x,y,z and 



ij 

fij = I ^ij I ) Iljj = , Vij = Tj Yj . 

fij 

The dipolar tensor enjoys the properties 

of which the first is exact and the second is asymptotically exact for a macroscopic sample 
with a large number of spins N ^ 1. The total magnetic field 

B = Boe, + (5i + H)e^ (6) 
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consists of an external longitudinal magnetic field Bq, transverse magnetic field Bi, and a 
feedback field H of the resonant electric circuit. In what follows, the longitudinal magnetic 
field is assumed to be constant and directed so that 



l^oBo < . (7) 

In particular, for electronic spins, Bq > 0, since /lo is negative. 

For what follows, it is convenient to pass to raising and lowering operators 

which are Hermitian conjugated with each other. Also, introduce the notation 

^\{ct^ + 2tC!f - Cf) . (8) 
Then the individual term (2) can be written as 

Hi = -iioBoS^ - D{Stf - \ ii^{B, + H) (st + S-) . (9) 
And the interaction term (3) takes the form 

Hij = ttij (^S-Sj - - S^S~^ + 2b*jSlSj + 2bijS~Sj + 

+ c^jStS+ + CijS-S- . (10) 
Thus, the spin system is described by the Hamiltonian (1), with the terms (9) and (10). 

5 Electric Circuit 

The external fields Bq and Bi are assumed to be given. What is not yet defined is the 
resonator feedback field H, produced by a resonant electric circuit. Let the circuit be char- 
acterized by resistance R, inductance L, and capacity C. The spin sample is inserted into 
a coil of n terns, length I and cross-section area Ac- The electric current in the circuit is 
determined by the Kirchhoff equation 

in which E is an electromotive force and $ is a magnetic flux 

$ = — nA^riM^ (12) 
c 
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formed by the x-component of the magnetization density 

M^ = ^Y.<St>, (13) 

with the brackets < . . . > implying statistical averaging. The filling factor rj is approximately 
equal to 77 ~ V/Vc, where V is the sample volume and Vc = AJ is the coil volume. 
The current, circulating over the coil, produces a magnetic field 

H = — J, (14) 

with c being the light velocity. Hence, Eq. (11) may be rewritten for the field (14). Introduce 
the circuit natural frequency 

Vlc \ cH j ^ ^ 

and the circuit ringing time 
with the related circuit damping 

where Q is the quality factor. Define the reduced electromotive force 



Then the Kirchhoff equation (11) can be transformed to the equation 

^ + 27// + f H{t') dt' = 7/1 - Attt) ^ (19) 
at Jo dt 

for the feedback magnetic field produced by the coil. 

The feedback equation (19) can be presented in another equivalent form [15-18] which is 
very useful for solving the evolution equations. For this purpose, we apply to Eq. (19) the 
method of Laplace transforms and use the transfer function 

G{t) = ^coscDt — ^ sincDt^ e"'''* , 
with CO = \/iJ^ — 7^. Then from the feedback equation (19) we obtain 

C Git - t') hhit') - AT^rjMJt')] dt' , (20) 

Jo L -I 

where = dM^/dt. 
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6 Evolution Equations 



First, we write the Heisenberg equations 



ih 



dSl 
dt 



St,H 



for the spin operators, using the commutation relations 

[Sr, Sj] = -25,^St , [St, Sf] = ±5,^St , [^r, [S'.f] = 5,, (S-St + StSr) . 
For the lowering operator, this gives 

dsr 



ih 



dt 



+ E 

mi) 



(21) 



The equation for the raising operator is obtained from Eq. (21) by Hermitian conjugation. 
The equation for S? is 

J CZ 1 



dt 



(22) 



To analyze these equations, we employ the scale separation approach [14,17,18,71]. No- 
tice, first, that in Eqs. (21) and (22), we may separate the terms 



I ^ /I 



h I2 "-^^ 



(23) 



which play the role of local fields acting on spins. Statistical averages of these fields, with 
the usage of the uniform approximation, are zero owing to the equalities 



E = E 







following from the properties (5). At the same time, statistical averages of the spin operators 
are certainly nonzero. Therefore the local fields (23), acting on a short scale, can be treated 
as operator variables of nature different from the spin operators. Such local fields may be 
modelled by random variables [68,69,72]. Thus, wc have two types of variables in the system, 
spin operators S^^ , S^, S^, and stochastic fields ,^oi ^7 The former are responsible for 
long-range global phenomena while the latter, for short-range local effects. The stochastic 
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fields describe fast fluctuations in the local surrounding of each spin. The existence of such 
fluctuations yields inhomogeneous dynamic broadening. 

To make the problem closed, it is necessary to define the stochastic averages over the 
random fields (23). The latter can be treated as white noise with the stochastic averages 

< Co{t) » = < m > = , < Ut)Ut') > = 2735(t - t') , 

« umt') » = « mat') » = o , « cmif) » = 2^sSit - to , (24) 

where 73 is the width of dynamic broadening. The method of modelling local fields by 
random variables can be called randomization of local fields. 

Averaging over spin operators, which are responsible for long-range phenomena, we may 
employ the mean-field approximation 

<srs;;> = <sr><sf> (z^j). (25) 

This kind of approximation can be used only for the pairs of spins at different sites i 7^ j, 
since and Sf commute in such a case. But the spin operators do not commute for i — j. 
Hence the averages < S°'Si > cannot be factorized as above. Notice that for spin S — 1/2, 
the anisotropy term in the Hamiltonian does not contribute to the equations of motion, 
since {Sf)^ = 1/4. The problem of spin decoupling, more general than condition (25), has 
been considered by several authors [73-76]. The term in Eq. (21), caused by the magnetic 
anisotropy, can be decoupled as follows: 

<SrSt + StSr>^l2- ^^<S-><St> . (26) 

This presentation enjoys correct limiting properties. For S = 1/2, it nullifies; while for 
5" — > 00, when spins behave classically, a simple factorization occurs. The decouphngs (25) 
and (26) do not take account of spin-spin correlations, which can be incorporated into the 
evolution equations by including the term describing spin attenuation, with the related spin- 
spin relaxation width 72. To allow for the influence of lattice, account must be taken of the 
spin-lattice relaxation, with the corresponding relaxation parameter 71. 

Let us average the equations of motion (21) and (22) over the spin operators, not touching 
the stochastic fields (23). In so doing, we introduce the following definitions. The variable 

x^^<Sr> (27) 

describes the rotation of transverse spin components. As is discussed in Section 3, this aver- 
age is connected with the arising transition coherence, or transverse coherence, or radiation 
coherence. The degree of such a coherence can be characterized by the real function 

y^j^ <St><Sr> = \x\\ (28) 

And the longitudinal spin polarization is given by 

z = ^ <St> . (29) 
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Recall that the radiation of spins happens at radio-frequencies whose wavelengths are much 
larger than the mean distance between spins. Therefore, it is admissible to employ the 
uniform approximation, assuming that the functions (27) to (29) do not depend on spatial 
variables. 

Note that, instead of resorting to the uniform approximation for the functions (27) to (29), 
it would be possible to work with the following arithmetical averages: transition function 

1 ^ 

coherence intensity 

1 ^ 

and longitudinal polarization, or spin polarization 

E < -5 • > . 



The resulting equations would be absolutely the same. 
Le us define the Zeeman frequency 

ujQ = ^1/^0^0 1 , (30) 

the effective transition frequency 

uo = uJo- {2S - l)Dz , (31) 

and introduce the notation 

/ = -^//o(Si + i/) + ^ (32) 

for an effective force acting on spins. Then from Eqs. (21) and (22), we obtain the evolution 
equations for the functional variables (27) to (29): 

dx 

— = -i (a;o + 6 - il2) x + fz , (33) 

^ = -272y {x*f + fx) z , (34) 

f^^-l{x*f + rx)-M^-a). (35) 
These equations are assumed to be comphmented by the initial conditions 

x{0) = xo , 2/(0) = yo , z{0) = zq . 
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Equations (33) to (35) are stochastic differential equations, since they contain the stochastic 
variables and ^. Getting stochastic equations is the price for making them closed. These 
also are nonlinear equations because of the resonator feedback field entering through the 
effective force (32). This resonator field is given by Eq. (20), where one has to substitute 
the magnetization density 

^^ = \ Pl^oS{x* + x) (^p=^^ . 

Due to the integral form (20), Eqs. (33) to (35) plus (20) compose a system of stochastic 
nonlinear integro-differential equations. 



7 Stochastic Averaging 

The system of equations (33) to (35) plus (20) looks rather complicated. Nevertheless, it can 
be essentially simplified by invoking the method of stochastic averaging [17,18,71], which is 
a generalization of multiscale averaging techniques to stochastic differential equations. The 
applicability of the method is due to the existence of several small parameters. 

The spin-lattice and spin-spin relaxation parameters are assumed to be small as compared 
to the Zeeman frequency, 

^ < 1 , ^ < 1 . (36) 
The interaction energy of spins with the resonator field is proportional to the natural width 

70 = ^ ■np^^ls . (37) 

Since 72 ~ UopplS"^ /h, where Uq is the number of nearest neighbors, then 

7o a 72 < 72 ■ 

nob 

Hence, the natural width (37) is small, as well as the width of dynamic broadening, 

^ < 1 , ^ < 1 . (38) 

Here we do not explicitly consider hyperfine interactions, surmising that their influence 
can be included in the values of the related relaxation parameters. In order to estimate 
the impact of these interactions, e.g. on the magnitude of the dynamic broadening, it is 
necessary to discriminate the cases when the radiating objects are nuclear spins or electronic 
spins. Recall that the spins of magnetic molecules are of electronic nature. The dynamic 
broadening of electronic spins, caused by dipolar and hyperfine interactions, respectively, is 
73 ~ PeA*e ^^cl 73 ~ pn^J'nPe, where pe and pn are the densities of electronic or nuclear spins, 
with pe = gsf^'sS and pn = giPNl being the electronic and nuclear moments. Similarly, 
denoting by capital letters the relaxation parameters for nuclear spins, we have the dynamic 
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broadening, due to dipolar or hyperfine interactions as Fa ~ PnA*n Fg ~ peHel^n- From 
here, the following relations are vahd: 

73 PnfJ'n pe/ie 



73 Pe/^e J- 3 Pn/^r 

Taking into account the equalities 

/^e _ 9sI^bS fiB _mp _ 



1836 , 



we see that the hyperfine interactions do not play an important role for electronic spins but 
can be rather important in the case of nuclear spins [28-30]. 
The external transverse magnetic field is taken in the form 

Bi = ho + hi cos ujt . (39) 

And let the resonant part of the reduced electromotive force (18) be presented as 

h = he cos ujt . (40) 

Introduce the notation 

''^^^ ^^ = ^' ^^=^- ^''^ 
The amplitudes of the fields (39) and (40) are supposed to be small, in the sense that 

^«1, ^«1, ^«1. (42) 

Uq Uq UJq 

The influence of stochastic fields has to be considered as weak, since the stochastic 
averages (24) are proportional to the dynamic broadening width that, according to Eq. (38), 
is small. In this way, the effective force (32) can be treated as weak. 

Finally, the resonant circuit is assumed to be of good quality, implying that 

^ < 1 (g > 1) . (43) 

The existence of the listed small parameters shows that the transition function x, defined 
by Eq. (33), is to be considered as fast, compared to the slow functions y and z, satisfying 
Eqs. (34) and (35). Conversely, y and z are temporal quasi-invariants with respect to x. 

The resonator field H , in the first approximation, may be found by iterating Eq. (20) 
with the solution of Eq. (33) of zero order with respect to small parameters, that is with 
X ~ Xo exp(— ia)oi), where 2; is a quasi- invariant. This iteration yields 

^ ^i{ax- a*x*) + 2/3cosa;i , (44) 
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where a is the couphng function of spins with the resonator feedback field, 



a = 7oa;o 



1 — exp{— — ujo)^ — 7O _|_ 1 — exp{i(ct; + ujo)t — jt} 



7 + i(a;-a;o) ^ — {{u + Uo) 

and (3 is the couphng function of spins with the electromotive force, 



(45) 



(46) 



Clearly, the action of the feedback field can be efficient only in the case of a resonant coupling, 
which requires the resonance condition 



|A| 

CO 



< 1 , A = a; - |a;o| 



(47) 



Note that the effective frequency (31) can, in general, be positive as well as negative. The 
spin-feedback coupling (45) simplifies if the resonance is good, which means that |A| < 7. 
In such a case, Eq. (45) acquires the simple form 



a 



770^0 
72 + A2 



(1 



-7* 



(48) 



More generally, the coupling function (45) is complex, with its real and imaginary parts 
being 



Re a — 



Im a — 



770^0 
72 + A2 

77o|'^o| 



cos At 



A 

7 



sin At e' 



-It 



A / ~ A ~ \ , 
(sin At + — cos An e~^* 



72 + A2 [7 

where the resonance condition (47) is used. 

Expression (44) is to be substituted into Eqs. (33) to (35). With this in mind, we define 

the collective frequency 

Q, = Cjo- zlm a (49) 

and the collective attenuation 

r = 72 - zRe q; . (50) 

These quantities depend on time through the slow variables z and a, because of which one 
may say that the frequency and attenuation experience dynamic shift. In the case of good 
resonance, when a is given by Eq. (48), expressions (49) and (50) are 



= ujq = ujq — {2S — l)Dz , r = 72 — q;^; . 
Also, we define the effective force 

/i = —iuq — 2i(i/i + (3) cos ut + ^ . 



(51) 



(52) 
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Then Eqs. (33) to (35) rearrange to 



dy 
dt 



dx 

— = + Co ~ '^r)a: + fiz — azx* , 

(JjV 



-2Ty + {x*fi + f^x) z-az {x*f + x 



^ = -ay - ^ {x*fi + fix) - 7i(z - a) + ^ a [{x*f + 



(53) 

(54) 
(55) 



The time derivatives of a and /3 are proportional to 7, because of which these functions are 
to be treated as slow, compared to x. Hence, a and /3, as well as Q and F, are temporal 
quasi- invariants with respect to x. 

The solution to Eq. (53), with the quasi-invariants kept fixed, reads 



+ z 



X = a;oe-('^+^)* exp |-i ^* ^0(^0 dt'^ + 
j^*/i(Oe"^'''^^^^*"*'^exp|-i^*eo(i") ti^"} ti^' ■ 



(56) 



The counterrotating term of Eq. (53), containing x*, gives to the form (56) a small addition 
of order 72/ci;o, because of which the latter is omitted. 
Define the quantities 

f = r + 73, 5 = uj-\n\, (57) 



which in the case of good resonance are 

r fa 72 — + 73 , 5 ~ a; — |a)o| 
Averaging Eq. (56) over the stochastic variable ^0, we get 



A . 



Xq + 



{1^1+ P)z 

5 + ir 



-{in+r)t 



The fast solution (56) has to be substituted into Eqs. (54) and (55) for the slow variables, 
whose right-hand sides are to be averaged over stochastic fields and over time in the infinite 
interval. Take the initial condition for the transition function x in the real form 

xo = < S^{0) > , 

which is not principal but just slightly simplifies the following formulas. Also, keep in mind 
that the collective attenuation is small, compared to the collective frequency. 



r 



< 1 
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From here, since 73 -C tuo, it follows that |r| <C And let us introduce the effective 
attenuation 

..,3-4.-^^^<^e-.<^(l-e-). (58) 

In the latter, the terms related to the action of the constant transverse field are of second 
and third order in small parameters. Omitting these terms gives 

Averaging the right-hand sides of Eqs. (54) and (55) over the stochastic fields and time, we 
come to the evolution equations 

^^-2{^^-az)y + 2Jz\ (59) 
dz 

— = -ay - Jz- 71(2; - a) , (60) 

whose solutions are called the guiding centers. 

As is evident, the derived equations (59) and (60) for the guiding centers are incompa- 
rably easier to analyse than the initial equations (33) to (35). Since Eqs. (59) and (60) 
are nonlinear, their solutions can display rather nontrivial behaviour, especially when the 
resonant transverse field or electromotive force are present. Spin dynamics, in the presence 
of such transverse injected fields is known to be quite complicated [13,77,78]. But even when 
external fields are absent, the nonlinear spin dynamics may demonstrate very interesting 
effects. 

Before analysing the evolution equations (59) and (60), let us pay attention to the physical 
conditions under which these equations have been obtained. The principal point here is the 
resonance condition (47). It is only when the spin system is in resonance with the coupled 
electric circuit that we could expect the appearance of noticeable transition coherence, which 
would develop in a self-organized way. This coherence, of course, can be induced by external 
transverse fields. However recall that our major aim is to consider superradiant regimes 
when the transition coherence and, respectively, radiation coherence, arises spontaneously, 
without being stimulated by intensive external fields. 

There can be several possibilities for reahzing the resonance condition (47). First of 
all, for systems with spin one-half, there is no magnetic anisotropy so that the effective 
transition frequency (31) coincides with the Zeeman frequency, cDq = cuq. Then we have the 
same situation as for nuclear magnets [14-20] with S = 1/2. 

If we are dealing with higher spins, then, nevertheless, there is the possibility of reducing 
the problem to an effective spin-one-half system. This can be done by tuning the resonant 
electric circuit to one of the transition frequencies of admissible 25" transitions and by sup- 
porting the population of the upper level with the help of a permanent nonresonant pumping, 
as it was accomplished for the nuclear spin / = 5/2 of ^''Al in experiments [11-13]. In that 
case solely the regime of pulsing spin superradiance can be achieved. 
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For higher spins, the transition from the state with the spin projection S to that with the 
projection —S corresponds to a multiphoton transition through 25" — 1 intermediate levels. 
If the longitudinal external magnetic field is sufficiently strong, such that (25" — 1)D <^ loq, 
then again cZ^o ~ ^O) ^^^d the considered 2S'-photon transition is not much different from the 
one-photon transition, with the same theory [14-20] being applicable. 

The worst case occurs when we are interested in a multiphoton transition in a spin sample 
with high magnetic anisotropy, such that (25" — 1)D » ujq. Then the effective transition 
frequency (31) is Oq ~ —(25' — l)Dz, which changes with time together with z. Because 
of this, to keep the resonance condition (47) valid, one should respectively vary with time 
the natural resonator frequency uj. It is feasible, in principle, to imagine such a sliding 
resonance, when the circuit characteristics, as inductance and capacity, are changing in time 
so that to preserve the approximate equality a; ~ (Dq, following the varying uq. But it looks 
that such a sliding resonance would be difficult to realize experimentally. However, even 
not this is the major obstacle. When the magnetic anisotropy is the prevailing part in the 
effective frequency (31), so that ujq ~ —(25* — l)Dz, then the spin-resonator coupling (48) 
is proportional to —Dz. As a result, az ~ —Dz^, which is always negative for the easy-axis 
anisotropy, with D > 0. Hence, the collective attenuation F, given in Eq. (51), is always 
positive. This means, according to Eqs. (53), (54), or (59), that the function y decreases with 
time. Therefore, there is no generation of coherent radiation. For this to occur, the collective 
frequency F must, at least for some period of time, become negative. As is evident from Eqs. 
(54) and (59), a negative attenuation F leads to the generation of coherent radiation, but a 
positive F leads to the decay of transition coherence. Thus, a too strong easy-axis magnetic 
anisotropy suppresses transverse coherence, hindering the development of multiphoton spin 
superradiance. 

In the case of an easy-plane anisotropy, with D < 0, the quantity az ~ |-D|z^ is positive. 
Hence, F = 72 — q:z could be negative. However, a strong x — y-plane anisotropy implies that 
it is difficult to noticeably polarize spins in the 2;-direction. That is, z is always small, which 
prevents making az large, such that > 72. Therefore a strong single-site anisotropy is an 
obstacle for achieving a well-developed transition coherence. 

8 Quantum Stage 

Let us assume that there are no strong external fields imposing on the spin system an 
essential transverse coherence, so that at the initial stage the motion of transverse spins is 
not coherent. This is the most interesting case to consider how the transverse coherence 
develops in a self-organized way from the initially incoherent motion. The period of time, 
when there are yet no collective correlations but only quantum spin interactions are present 
can be called the quantum stage. 

At the very beginning of the process, when ■jt ^ 1, the spin-resonator coupling function 
(48) is yet close to zero, and we may set a — > 0. If there are no transverse external resonator 
fields, which imphes that ui — Uc — 0, hence P — 0, then the effective attenuation (58) 
reduces to J = 73. In that case, when collective effects, due to the correlation of spins 
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through the feedback field have not yet been developed, the evolution equations (59) and 
(60) are yet quite simple, having the form 

dy ^ ^ o dz . . 

— = -272^ + 273^ , = -(71 + 73)^ + 710- ■ (61) 

Then at short time t ^ 0, we have 

y^L- + z^L- ^) + ^ . (62) 

V 72 y 72 V 71 + 73 / 71+73 

With time, spin correlations increase owing to the resonator feedback field. Strengthening 
spin correlations result in the developing transverse coherence. When the transition coher- 
ence is well developed, collective phenomena come into play and the dynamic behaviour 
of spins becomes qualitatively different from that at the initial stage. The change in the 
features of spin motion occurs, of course, gradually. However, it is possible to define the 
moment of time separating these two regimes of motion. The qualitative change in spin 
dynamics occurs when the collective attenuation width (50) changes its sign because of the 
growing spin-resonator coupling. Then the quantum stage transfers to the coherent stage. 
The crossover time tc, separating these stages is given by the equation 

TiQ = . (63) 

The latter, in the case of good resonance, takes the form 

«(tc) z{Q = 72 • (64) 

Introducing the effective spin-resonator coupling parameter 

^^nou^ (65) 
^ 72(7' + A2) ^ ' 

in which 

ujq = cuo — {2S — 1)Dzq , A = a; — pol , (66) 
we obtain from Eq. (64) the crossover time 

t.^rln(^]. (67) 

The latter, to be positive and finite, requires that 

gzo > 1 (0 < tc < 00) . (68) 

When condition (68) is not valid, then the quantum stage will be never replaced by the 
coherent one. For instance, if gzo — 1, then tc ^ 00. When the initial spin polarization is 
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sufficiently high and the spin-resonator couphng (65) is strong, then the crossover time (67) 
is 

tc- — {gzo > 1) . 

For the coherent regime to appear, it is necessary that the crossover time (67) be shorter 
than the following relaxation times, 

t,^lT, = -,T, = -,n = -] . (69) 
[ 7i 72 73 J 

In the other case, spin polarization would relax before the transition coherence could arise. 
Then the solutions (62) at the crossover time can be written as 

y{tc) ^ ?/o + 273tc^o > ^(Q '^zo + jitca . (70) 

If the single-site anisotropy is strong, such that (25* — l)Dzo ^ luq, then ujq ~ —zq. Since 
the spin-resonator coupling (65) is proportional to (Dq, then g ~ —zq, from where it follows 
that gzo ~ —Zq, which is always nonpositive. Hence, condition (68) cannot be hod true. 
This means that the coherent stage can never appear. In order that the transition coherence 
could develop, the external magnetic field Bq has to be sufficiently strong so that to suppress 
the destructive role of the magnetic anisotropy. 

The dynamic broadening width 73 in Eqs. (61) has been assumed to be a constant. In 
general, it can be a function of time. For example, considering the dynamics of the spin 
polarization at the beginning of the process, when t ^ 0, and when 71 <^ 73, we have the 
equation 

dz , 

^ = -T3. (t-0). 

Setting here 

^'^'^ " 2(i^ ^ 
results in the square-root relaxation [66] 

Z ~ Zo{l - y/^) 

observed in molecular magnets [45] at short times, below the blocking temperature. Thus, 
for Mni2 one has 7s ^ 10~^ s. 

At longer times, below the blocking temperature, the relaxation of magnetization in 
molecular magnets follows a stretched exponential law [37,42,44,49,51,52]. This can be de- 
rived if 

73(t) = ^ (7*^)'' (0 < «: < 1) , 

which gives 

z = 2;ocxp{-(7*t)'^} . 

The power k ^ 0.5 for T < Tb and k = 1 for T > T^. The relaxation parameter 7* depends 
on the applied magnetic field. 
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Recall that all the consideration of this section concerns the situation when collective 
effects have not yet come into play, so that the spin-resonator coupling is negligible, o; ~ 0. 
But if the spin system is coupled to a resonator and condition (68) is valid, then there exists 
a finite crossover time (67) when the quantum stage is followed by the coherent one. 



9 Coherent Stage 

After the crossover time (67), the coupling function (48) fastly increases to the value a ^ 
g'f2- This means that collective effects come into play and the transition coherence starts 
developing. At the transient stage, when time is larger than tc but much shorter than Ti, one 
may neglect the longitudinal relaxation, omitting 71 in Eq. (60). Superradiance occurring 
at this transient stage, when 

tc<t<^Ti, (71) 

is the transient superradiance. 

If there are no transverse external fields, then the effective attenuation (58) is J = 73. 
For a sufficiently large coupling parameter (65), such that (772 ^ 73, one may also neglect 
in Eqs. (59) and (60) the terms containing 73. Then these equations read 

dy dz 

— = -272(1 - gz)y , — = -^72y . (72) 



y = ( ^Vsech^ ( ^—^] , z^- ^tanh ( ^ I + - , (73) 




Equations (72) can be solved exactly [15-18] yielding 

sech^ ( - — — ] , z — — 

912 V 

where the pulse width 7^, the pulse time r^, with the relation '^pTp = 1, and the delay time to 
are the integration constants to be defined from the initial conditions taken at the crossover 
time tc- Equating the functions (73) to the values (70) gives for the pulse width the equation 

7p = 7g + {9l2f{yQ + 273*c^o) ' Ig = 72(1 - qzq) , ^pTp = 1 , (74) 
and for the delay time, we get 

7p-79 



to^tc+^ln 



Ip + lg 



(75) 



Solutions (73) describe a superradiant pulse. Since, by definition, superradiance is a 
self-organized process, we consider the case when there is no strong transverse coherence 
imposed on the system at the initial time. This implies the inequality 

^'yo < 1 . (76) 

In order that the delay time (75) would be finite, and larger than the crossover time, so that 

< to < 00 , (77) 
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it is necessary that 

yo + 2-f3tczl > , gzo>l. (78) 

One more condition on superradiance is that the pulse time be shorter than the spin-spin 
dephasing time, 

Tp < T2 . (79) 

Then, taking into account that 73^0 <S 1, we may distinguish two superradiant regimes: 
triggered superradiance, when 



gzo>l + ^Jl-g^yo, Vo^O, (80) 
and pure superradiance, when 

gzo>2, 1/0 = 0. (81) 

One of the necessary conditions for superradiance to occur is gzo > 1, which is the same 
condition (68) for the existence of the crossover time (67). This condition, in view of Eqs. 
(65) and (66), can be written as 



uozo > {2S - l)Dz^ + 



2 , 72(7' + A') 



77o 



which tells us that the longitudinal magnetic field Bq has to be sufficiently strong and the 
initial spin polarization must be positive, Zq > 0. 

For the pure spin superradiance, when yo — 0, taking into consideration that ^y^tc <^ 1, 
the pulse width and time can be presented as 

T 

7p = {gzo - 1)72(1 + 73^c) , Tp = ^ (1 - 73tc) , (82) 

gzo - 1 



and the delay time (75) takes the form 

to^tc+^ In 

where the notation 



73^c 



(83) 



U^Ue'-^^^^rl^^Un ( 

gzo-l \gzo-l 



is used. When <^ 1, which is equivalent to the inequality gz^ ^ 1, then ic — tc- 

At the delay time t — to, the coherence intensity y is maximal, which follows from the 
solutions (73) giving 

y{to) = (^zo - ^ j (1 + 273^,) , z{to) = ^ . (84) 
After this, the transition coherence fastly decays. Thus, for t ^ to, one has 

y ~ Ay{to)e-^^-' , z ~ -zq + ^ . (85) 
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For large g ^ 1, the initial spin polarization becomes almost completely inverted. This 
effect can be used for fast polarization reversal in polarized solid-state targets employed 
in scattering experiments [10,79]. The polarization reversal under a superradiant burst is 
illustrated in Fig. 1, which results from numerical simulations [21-23]. 



10 Pulsing Regime 

After the transient superradiant burst, occurring at the delay time to, the transition coher- 
ence dies out, unless the inversion of spin polarization is supported by a permanent pumping. 
This pumping can be accomplished by means of dynamic nuclear polarization. In the latter 
case, the regime of pulsing spin superradiance can be achieved. To study this regime, we 
consider a spin system without external transverse fields. Then the effective attenuation (58) 
is J = 73. At long times, such that ^ 1, the coupling function (48) reaches its maximal 
value a = g'j2- In that case, Eqs. (59) and (60) can be written as 

dy dz 

with the right-hand sides 

vi = -272(1 - gz)y + 273^^ , V2 = -gi2y - 1-iZ - 71(2; - a) . (87) 

In the presence of the permanent polarization pumping, the pumping rate 7^ takes the place 
of the spin-lattice relaxation parameter 71, while a G [—1, 1] is a pumping parameter. 

To understand the behaviour of the solutions y and z at long times t — 00, let us find 
the stationary points of Eqs. (86), which are given by the equations V1 — V2 — 0. In general, 
these equations give us two stationary points y^, z^ and 1/2, Z2, also called fixed points. 
Among these, however, only those have sense, which satisfy the physical restrictions 

0<|/<1, -l<z<l . (88) 

Also, the actual solutions of Eq. (86) will tend to that fixed point which is stable, whose 
properties can be derived by means of the Lyapunov stability analysis. The stability of fixed 
points depends on the values of the spin-resonator coupling (65) and of the effective pumping 
parameter 

a* = . (89) 

7! +73 

The value of this parameter depends on the relation between the pumping rate 7^" and the 
dynamic broadening width 73. Since the effects of temperature can be modelled through a 
thermal bath of random fluctuations [80] , the temperature dependence can be incorporated 
in 73. 

When the pumping is weak, such that 

ga* < -1 , (90) 
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then the stable fixed point corresponds to the solutions 

(7^ + 73)73 M , *| * , ^3 \ /^-.x 

Vi- r-* \^ + 9(^\, ^i~(T 1 + -] . 91 

The related characteristic exponents are 

X+ ~ -7* - ll^lJl ^3 ^ ~ -272(1 - ga) - 73 . (92) 

9(^12 

The fixed point (91) is a stable node. The coherence intensity yl ~ 73/5'72 is small. 
When the spin-resonator coupling is weak, so that 

\ga*\ « 1 , (93) 

then the stationary point is given by the solutions 

* 7i73(g*)^ f, , 7i-73 *\ * 73 *\ i^.. 

Vi - T^, ^ 1 + ~r, ] , ^1 ~ (7 1 — ga . (94) 

(7i + 73)72 V 7r + 73 / V 7i + 73 / 

The corresponding characteristic exponents 

^1+ ^ -71* -73 (1 - . !^'\ 9'J*] , ^ -272 (1 - 1'^'''^' 9'^*] (95) 

V 7i - 272 + 73 / V 7i - 272 + 73 J 

show that the fixed point (94) is, as early, a stable node. The hmiting value <^1. In these 
both cases of either Eq. (90) or Eq. (93), the solutions monotonically tend to the stationary 
point y*, zl. 

When the spin-resonator coupling as well as pumping are sufficiently strong, so that 

ga*>\, (96) 
then the solutions of the evolution equations (86) tend to another fixed point 

1/2*^^^(5^*-!), ;.;^-fl-^). (97) 

r72 9 \ 9<^ii) 

Now, the characteristic exponents are complex valued. 



Xf^-i7,--i73(l + ^)±i^», (98) 



containing the asymptotic frequency 



2gaYil2 . (99) 



The fixed point (97) is a stable focus. The limiting value of the coherence intensity 1/2 ~ 
I1/9I2 is small. However, in the way to the fixed point (97), the spin system exhibits a 
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series of superradiant bursts, which can be demonstrated by solving Eqs. (86) numerically 
[20,30,71]. This is the regime of pulsing spin superradiance. The superradiant bursts are not 
equidistant in time, and become approximately periodic only at t — > oo, with the asymptotic 
period 



where 7*Tj* = 1. Each superradiant burst, occurring in the time interval < t < T^, 
possesses a high level of transition coherence, but after the amplitude of the coherence 
intensity y diminishes, tending to 1/2. 

In conclusion, superradiant regimes may be employed in the operation of spin masers. 
There are three possible such regimes: triggered superradiance, pure superradiance, and 
pulsing superradiance. In materials with higher spins, as molecular magnets, it is necessary 
to take into account single-site magnetic anisotropy. If this anisotropy is strong, it hinders 
the possibility of achieving a high level of transition coherence in multiphoton transitions. 
However, by a sufficiently strong external magnetic field the destructive role of magnetic 
anisotropy can be suppressed. Another possibility is to tune the resonant electric circuit to 
one of the quantum transition lines of a high spin. When a sole transition line is selected 
by the resonant tuning, the situation becomes similar to the case of an effective two-level 
system 
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Figure Caption 

Fig. 1. Intensity of transition coherence y in arbitrary units (upper curve) and the 
longitudinal spin polarization z (lower curve) as functions of time measured in units of T2, 
obtained from numerical simulations for 300 spins. 
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